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ABSTRACT 

We propose here a robust scheme to infer the physical parameters of compact stars from their f- 
mode gravitational wave signals. We first show that the frequency and the damping rate of f-mode 
oscillation of compact stars can be expressed in terms of universal functions of stellar mass and moment 
of inertia. By employing the universality in the f-mode one can then infer accurate values of the mass, 
the moment of inertia and the radius of a compact star. In addition, we demonstrate that our new 
scheme works well for both realistic neutron stars and quark stars, and hence provides a unifying way 
to infer the physical parameters of compact stars. 

Subject headings: gravitational waves — stars: neutron — stars: oscillations (including pulsations) — 
equation of state — relativity 
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1. INTRODUCTION 

Ever since the fi rst pr ediction of neutron stars (NS) by 
iBaade fc Zwickyl |l934), NS have become a major test 
bed of various theories of dense matter and nuclear mat- 
ter. The extreme high pressure and density achievable 
inside NS offers an ideal and unique environment to infer 
the equation of state (EOS) of nuclear matter from as- 
tronomical observations. On the other hand, the possible 
existence of quark stars (QS) would provide a strong sup- 
port to a new form of ma t ter, namely, t he quark matter 
(|Itohlll970L rBodmeilfl97ll I Wittenl fl98l . Moreover, the 
size and the mass of QS could also cast light on the prop- 
erties of quark matter. Thus, extracting useful informa- 
tion and constraints about nuclear and particle physics 
from astronomical observations of compact stars like NS 
and QS has become an active direction in astrophysics 
(see, e.g.. lLattimer fe Prakashll2007l lHaensel et alJl2007l 
for reviews). 

Since the density in the core of a NS (or QS) can be 
several times the normal nuclear density, the EOS there 
is rather uncertain. Various physical phenomena includ- 
ing pion condensation, hyperon formation and decon- 
fined quark matter have been considered by theoretical 
nuclear physicists and sophisticated techniques of many- 
body calculations have been developed to yield a bunch 
of EOS intended for the description of the deep interior of 
compact stars. It is standard practice for astrophysicists 
to construct compact stars with these EOS for nuclear 
matter and establish various relationships among differ- 
ent physical quantities which could be obtained from as- 
tronomical observations. Motivated by the possibilities 
of constraining EOS of nuclear matter and extracting 
physical parameters of compact stars as well, researchers 
in astrophysics and nuclear physics have been ac- 
tively seeking for EOS-dependent and EOS-independent 
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relationships by examining the physical characteris- 
tics of co mpact stars constructed with different EOS 
(see, e. g.. lLattimer fc Prakashl 120011 IBeieer fc Haensell 



20021 lLattimer fc Schutzj |2005l iBeiger et all 120051 
Lattimer fc Prakashl I2007L lHaensel et alJl2007T K For ex 



ample. lLattimer fe Prakashl (|200ll ) compared the struc- 
ture of compact stars of various kinds and discovered 
several empirical relationships connecting different phys- 
ical characteristics of a star, such as the mass, the radius, 
the moment of inertia, and the mass distribution func- 
tion as well. Such relationships can be applied to infer 
the physical attributes of a compact star, including its 
EOS, from astronomical observations. 

The main objective of the present paper is to unveil 
the universality embedded in the pulsation frequencies 
of non-rotating compact stars and to exploit these find- 
ings to infer the physical characteristics of NS or QS (e.g., 
mass, radius, moment of inertia and EOS) from their os- 
cillation spectra. The pulsations of compact stars are an 
interesting and popular topic in its own right, for gravi- 
tati onal waves (GW) can be genera ted in such processes 
(see[Th orne fc Campolattaroll 19671 for a seminal exposi- 
tion of the topic). It is indisputable that the detection 
of GW would be a major milestone in general relativ- 
ity and astrophysics. Several Earth-based GW interfer- 
ometric detectors such as LIGO, VIRGO, GEO600 and 
TAMA300 have been operating. While the current detec- 
tors are still not sensitive enough to detect GW directly, 
interesting upper limits have been placed on either the 
GW strains or event rat e for several potentia l astrophys- 
ical sources (see, e.g., lAbbott et al.l l2009aHbl l3. for the 
results obtained from the latest science runs of LIGO). 

Owing to the energy carried away by GW, pulsations of 
compact stars are damped harmonic oscillations, which 
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Each QNM is character- 
ized by a complex eigenfrequency u) = uj r + ioji and 
has a time dependence exp(zwt), displaying exponen- 
tial decay with a damping time r = l/wj. Several 
attempts have been made to relate w r and oj; (or r) 
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to the mass M and the radius R of a compact star 
and some universal relationships which are, to certain 
degree of accuracy, EOS- independent have been found 
(lAndersson fc Kokkotaill996l. 119981 iBenhar et"al]|1999l 
12004 iTsui fc Leung 2005b). However, most of these re- 
lationships fail to describe QNM of QS, which are stiff 
and self-bound. 

On the other hand. lLattimer fe Schutzl (|2005h found an 
empirical relation approximately expressing the moment 
of inertia in terms of the mass and the radius. Such a 
relationship is universal as long as the EOS under consid- 
eration is not too soft (e.g., hyperon matter) or too stiff 
(e.g., quark matter). They proposed to use the univer- 
sality to estimate the radius of a compact star from its 
mass and moment of inertia, with the latter two quanti- 
ties being measurable for double pulsar systems by con- 
sidering the spin-orbit cou pling effect in general relativity 
(jLattimer fc Schutzll2005f ). 

Mot ivated by the discovery of lLattimer fe Schut3 
(2005), we propose here to investigate the correspon- 
dence between the QNM frequency (both ui x and a;;) of 
the (fundamental) f-mode oscillations, M and the mo- 
ment of inertia / of compact stars. We show that there 
exist EOS-independent universal relationships between 
them, namely equations ([6]) and (O, which are accu- 
rate up to a few percent or better and work nicely for 
QS as well. Based on equations ^ and ([7]), a feasible 
scheme is in turn established to determine the mass and 
the moment of inertia of a compact star once its f-mode 
frequency is found from GW observations. As the ra- 
dius R of a NS (or QS) is approximately expressible in 
terms of the mass and the moment of inertia, our scheme 
can also yield good estimate of t he radius and hence im- 
posing constraints on th e EOS (Bei ger fe Haense3l2002L 
lLattimer fc Schutd [20051 ) . 

The outline of the paper is as follows. Section 2 is a 
brief review on the the behavior of f-m ode oscillations of 
compact stars (|Tsui fc Leund [20 05b) and the universal 
relationships among the mass, the radius and the mo- 
ment of inertia of compac t stars (jBeiger fc Haensell2002L 
lLattimer fc Schutzl 120051 ). In Section 3 we introduce an 
accurate universal relation between the scaled frequency 
Muj and the physical quantity -JM 3 /I. In Section 4 we 
establish a feasible and robust scheme to apply our find- 
ings reported in Section 3 to invert the mass, the moment 
of inertia, and in turn the radius of a compact star if its 
f-mode GW signal is detected. We finally summarize our 
paper in Section 5 with a brief discussion. Unless oth- 
erwise noted, we adopt geometric units with G = c = 1 
throughout the whole paper, and use kilometers as the 
unit of lengths. 

2. GW AND MOMENT OF INERTIA OF COMPACT 

STARS 

As mentioned in Section 1, the oscillation modes 
of compact stars are damped due to the emission 
of GW. Hence, they are called quasi-normal modes, 
which can be further divided into polar and axial 
classes in each angular momentum sector according to 
the property of relevant metric and fluid perturba- 
tions under inversion. For typical non-rotating com- 
pact stars, the polar class includes (fundamental) f-mode, 
(gravity) g-mode, and (pressure) p-mode and (space- 
time) w-mode, while for non-rotating perfect fluid stars 



the axial class contains o nly the w- modes (see, e.g., 
iKokkotas fc Schnfi dt 1999, for a review on the prop- 
erty of these modes). The methods for evaluation of 
QNM of compact stars are well documented in several pa- 
pers and reviews (see, e.g., iLindblom fc Detweiler1ll983L 
Detweiler fc Lindbloml 119851 IC handrasckha r fc Ferraril 



1991. lKokkota7fc SchmidtJIW9^ 

The study of QNM of compact stars has a 
long history dating back to the pioneering work of 
iThorne fc Campolattarol (|1967l ) . It is by now well estab- 
lished that the study of GW mode spectra of compact 
stars could provide useful inform ation about the internal 
struc t ure of the stars (see, e . g.. lAndersson fc Kokkotai 
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that the direct detection of GW will be realized in the 
near future, different proposals have been put forward 
to determine the physical parameters (e.g., mass, radius 
and EOS) of a compact star from its GW spectra. In 
the present paper we focus our attention on an EOS- 
independent universal behavior of the quadrupolar (an- 
gular momentum index I = 2) f-mode and show that 
such universality can lead to accurate determination of 
the mass and the moment of inertia of a compact star 
(NS or QS). Combining our results with equations ©, 
([31) and ([SJ, one can also find an approximate value of 
the radius. 

The reason why we place special emphasis on the f- 
mode is that the relatively low frequency (~ 1 — 3 kHz) of 
the f-mode makes it more promising to be detectable by 
GW detectors than other modes. For example, consider - 
ing a compact star 10 kpc from us. IKokkotas et all (|200lD 
estimated that the energy required in the f-mode in order 
to lead to a detection with signal-to-noise ratio of 10 by 
the LIGO detector is 4.9 x 10~ 5 M Q c 2 . For comparison, 
in order to achieve the same signal-to-noise ratio, the en- 
ergy in the p-mode (w-mode) must be 4.0 x 10 _3 M Q c 2 
(6.8 x 10~ 2 Mqc 2 ). The energy requirement is much re- 
laxed for the advanced LIGO detector and the same 
signal-to-noise ratio is achieved if the energy in the f- 
mode is 8.7 x 10~ 7 M Q c 2 . Since the total energy released 
in violent astronomical events such as supernova explo- 
sions coul d be as much as 10~ 2 M ^c 2 in some optimistic 
estimates (K okkotas et al. 2001), it seems plausible that 
the f-mode QNM signal of a fiercely pulsating compact 
star inside our galaxy could be detected by future GW 
detectors, especially when the narrow-banding capability 
is available to enhance the sensitivity of the detectors in 
the high-frequency region at around 1 kHz. 

There have been several attempts to search for an ap- 
proximately EOS-independent behavior of the f-mode 
(lAndersson fc Kokkotas1ll996l 119981 IBenhar et^L1ll999l 
12004 ITsui fc Leundl2005bi rin particular. ITsui fc Leund 
( 2005bj) showed that the scaled frequency Mw for the f- 
mode of realistic NS is essentially a universal function of 
the compactness M/R, which can be approximated by 
the quadratic formulae: 



Mw r = 0.15 



Muj-, 



0.56 



0.02. 
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These two formulae are good approximation for most NS. 
However, the f-modes of QS do not follow equations (J]) 
and @. Instead, they are approximated better by an- 
other quadratic fit. The difference between the behaviors 
of the f-modes of NS and QS are expected because the 
density profile of a QS is nearly uniform, which is very 

different from that of a NS. 

On the other hand. IBeiger fc Haensell ((2002) and 
lLattimer fc Schutzl (|2005() discovered some empirical re- 
lationships relating I an d compactness M /R. For NS, 
IBeiger &; Haensell (|2002D found that 



8x10 



f _ I 1/(0.1 + 2a?) 
1 ~ \ 2(l + 5a:)/9 



if x < 0.1, 
if x > 0.1, 



(3) 



where the scaled mome nt of inertia I = I /MR 2 and 
x = (M/M Q )(km/R). lLattimer fc Schutzl (|2005D also 
proposed an alternative expression: 



I = 0.237(1 + 4.2a; + 90x 4 ) 



(4) 



which yields better agreement for NS with higher com- 
pactness. However, both equati ons (1H) and flU) fai l to de- 
scribe the case of QS. Instead, iBejger fc Haensell (|2002f ) 
found another relationship for QS: 



/=-(! + *) . 



(5) 



The distinction between the behaviors of NS and QS is 
also expected and attributable to the difference in the 
stiffness between nuclear and quark matters. 

Assu ming that both I and M are obtained from obser- 
vations, lLattimer &: Schut3 (|2005| ) further argued that 
the radius R can be found by solving these equation. 
However, two remarks about this method are in order. 
First, for the NS case the data points corresponding to 
stars characterized by different EOS are rather scattered 
and may deviate appreciably from equations (|3]) and ((4]) 
. As a consequence, the value of R obtained in this way 
is subjected to uncertainties depending on the EOS. Sec- 
ond, additional information has to be gathered to distin- 
guish between the data of QS and NS beforehand. Re- 
gardless of these remarks, it is observed that for a fixed 
EOS the scaled moment of inertia increases monotoni- 
cally with the compactness, indicating a one-one corre- 
spondence between these two quantities. In the following 
discussion we shall consider the feasibility to replace the 
compactness used in equations (JTJ) and @ with the mo- 
ment of inertia. 

3. UNIVERSALITY IN F-MODE OSCILLATIONS 

It is interesting to note from the reviews sketched 
in Sections 2 that both the scaled moment of inertia 
I /MR 2 and scaled f-mode frequency Mui are rather in- 
sensitive to EOS and are, withi n certain accuracy, uni- 
versal functions of compactness (Bcjgcr & Ha enselll2002L 
lLattimer fc Schutzll2005llTsui fc Leung||2005bh . Such co- 
incidence is interesting, reflecting the fact that compact- 
ness is a useful parameter characterizing a compact star. 
However, it is also noticed that the scaled f-mode fre- 
quencies (or the scaled moments of inertia) of NS and QS 
follow two different universal behaviors. The disparity in 
the universal behaviors of NS and QS is attributable to 
the difference in their density profiles. 
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Fig. 1. — The scaled frequency of f-mode QNM is plotted against 
the effective compactness r\. The solid lines in the lower and upper 
panels respectively represent the best quadratic fits in r\ to Mui r 
and Afa>j. In particular, the solid line in the lower panel is given 
by equation JB). 

In the present paper we propose to examine the possi- 
bility of using the moment of inertia as an independent 
parameter, instead of the compactness, to characterize 
the f-mode spectra and in turn to search for any uni- 
versality underlying these spectra. The physical motiva- 
tion of our proposal is that the moment of inertia carries 
richer information about the mass distribution inside a 
star, which certainly affects f-mode oscillations. Hence, 
the relationship between the moment of inertia and the f- 
mode is likely to be more direct and less EOS-dependent. 
Besides, as mentioned in Section 2, for compact stars 
with a fixed EOS, the moment of inertia bears a one-one 
correspondence with the compactness. To this end, we 
replace the compactness M/R by another dimensionless 
factor i] = y/ M 3 /I, which is proportional to the com- 
pactness for stars with uniform densities and is termed 
as the effective compactness. 

The scaled f-mode frequencies for various realistic NS 
and QS are plotted against the effective compactness rj 
in Figure [TJ where the f-mode QNM frequencies of the 
compact stars constructed with the following EOS are 
shown: 

1. Nine popular ordinary n uclear matter EOS fo r 
NS, includi ng Model A (IPandharipandd Il971al k 



Model C dPandharipandd 11971 



irio 
Lb): 



(|Canuto fc Chitrd 11974ft : Model I 



Model G 
(jCohen et all 
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TABLE 1 

Percentage error of different universal formulae 



quantity 


equation 


01 


0"2 


CT3 


OJ r 


m 


5.8 


4.2 


21 


OJ r 


m 


0.37 


0.56 


0.70 


Ul\ 


121 


220 


6.1 


55 


OJj 


w 


1.5 


2.1 


1.1 
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Fig. 2.— LOiI 2 /M 5 is plotted against q 2 (i.e., M 3 /I). The solid 
line is the best linear fit to the distribution (see equation H7} ). 



Il970h: APR 1 (lAkmalet alJflQQSh : AU, UU and 
UT jWiringa et al.ll988j) : and GM24 (iGlendennind 
11993, p. 244), which cover a wide range of softness 
(with an adiabatic index 7 ~ 1.6 — 4 in the high- 
density region). 

2. Two QS models, QS1 and QS2, constructed 
with the MIT bag model where the effects of 
strange quark mass and temperat ure are consid- 
ered (see, e.g.. IGlendenninel 11997( 1. Characteris- 
tic physical parameters for QS1 (QS2) are as fol- 
lows: the bag constant B is given by B 1 / 4 = 
154.5 MeV (150.0 MeV), strange quark mass m s = 
150 MeV (100 MeV), and temperature T = MeV 
(10 MeV). 

It is clearly seen that the data for NS and QS now 
display almost identical universal behaviors with high 
degree of accuracy. In fact, a quadratic expression in 77 
of the form 



MUJ V 



-0.0047 + 0.133? ? + 0.0575?7 2 



(6) 



suffices to yield a perfect fit to M w r of all stars (NS and 
QS) . On the other hand, despite the fact that Mu)\ of all 
stars fall on a common smooth curve, quadratic fits in rj 
do not work very well in this case. 

Motivated by the quadrupole radiation formula of 
GW and hence the sugges tion that uii oc M 3 /R 4 
(jAndersson fc Kokk otas 1998J) . we plot uj-J 2 /M 5 against 



rj 2 (i.e., M 3 /I) in Figure [5] and see that all data points 
now fall on a single straight line with high accuracy: 



7 2 Wi/M 5 



-0.00694 + 0.02567? 2 



(7) 



Therefore, we conclude that moment of inertia is a good 
physical parameter to unveil the universality inherent in 
the f-modes of compact stars. 

To compare the accuracy of various proposed universal 
curves for the f-modes quantitatively, we evaluate the 
root-mean-square percentage error a defined by 



^ N 



(8) 



Here N is the total number of f-mode data considered in 
the sample, yi corresponds to the exact value of the phys- 
ical quantity y (e.g., Mw, and Mu\) of the «-th f-mode 
data, iji is the approximated value of y inferred from the 
universal curves discussed above. In Table Q] we show 
the percentage error a in cj r and cjj for the data obtained 
from different universal curves and different samples of 
f-mode data, namely tho se of (i) NS shown in Figure [1] 
(<7i); (ii) NS considered in lAndersson fc Kokkotas! (| 19981 ) 
(172); and (iii) QS shown in Figure [T] (03) . 

It is clearly shown in Table [T] that the universal curves 
((6|) and ((7|) proposed in the present paper, which consid- 
ers mass and moment of inertia as two crucial parameters 
characterizing f-mode pulsations of compact stars, work 
nicely with high precisions. For w r , the error is less than 
1%, while for Wj, the error is only about 1-2%. In partic- 
ular, such high degree of accuracy prevails for both NS 
and QS, whereas the universal curves using conventional 
compactness fail to yield satisfactory result for QS. Thus, 
under the assumption that the f-mode GW signals from 
a compact star could be detected, we can then infer the 
mass and moment of inertia of the star from its f-mode 
using the universal curves ((6|) and ([7|), without any prior 
knowledge of its EOS and composition. We shall estab- 
lish a feasible inversion scheme to achieve such purpose 
in the following section. 

Before introducing the inversion scheme, we per- 
form an independent test here to show that equa- 
tions © and © are robust and accurate for other 
commonly considered EOS. We consider the f-mode 
QNM of NS constructed with five modern EOS o f 
nuclear matter, namely N 1H1 (iBalberg fc Gall [1997), 
BBB2 dBaldo et all Il997lh BPAL12 (iBombacil Il995l i. 
FPS (iPandharipande fc Ravenhalll U989D and SLy4 
(|Douchin fc HaensellteOOlE which have not been used in 
the derivation of equations ([6|) and J7|). The root-mean- 
square percentage errors a for these f-modes are just 
0.25% and 0.95%, respectively. This once again demon- 
strates the robustness of our finding reported here. 
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TABLE 2 

Percentage error of our inversion scheme 



EOS 


M 


SM/M 


6R/R 


81/ 1 


AU 


0.8 


-0.056 


-2.6 


-0.031 


AU 


1.0 


-0.13 


-0.67 


-0.16 


AU 


1.6 


4.0 


2.1 


6.9 


APRl 


0.8 


-0.066 


-5.1 


-0.23 


APRl 


1.2 


-0.21 


-1.6 


-0.19 


APRl 


1.6 


0.50 


-0.11 


0.78 


EOS A 


1.535 


0.40 


-1.2 


0.27 


EOS A 


1.328 


1.1 


-0.46 


1.8 


EOS B 


1.405 


-0.17 


-3.5 


0.54 


EOS B 


0.971 


1.4 


-4.1 


2.7 


GM24 


1.536 


-2.2 


-4.9 


-2.6 


GM24 


1.405 


2.5 


-1.1 


5.1 


QS1 


0.8 


-0.096 


18 (0.53) 


-0.18 


QS1 


1.4 


0.93 


12 (1.3) 


1.3 


QS2 


1.0 


-0.17 


17 (0.72) 


-0.3 


QS2 


1.6 


1.6 


11 (1.7) 


2.5 



4. INVERSION SCHEME 

Assuming that both uj r and lo\ of the f-mode of a com- 
pact star have been obtained from GW observations, we 
can then solve the simultaneous equations and ([7]) to 
determine the mass M and the moment of inertia / of 
the star. 

Moreover, making use of equation @ proposed by 
lLattimer fc Schutzl (|2005D . we can estimate the radius 
R with the input of M and /. Since all these three equa- 
tions hold only approximately, the inferred values of M, 
R and / are expected to deviate from the exact values by 
SM, 5R, and 51. As shown in Tabic [1] percentage errors 
in w r and u\ predicted by equations ([6]) and ([JJ are small 
and usually less than a few percent. It is reasonable to 
expect that the percentage errors 8M/M and SI /I are 
also small and such inversion scheme can lead to accurate 
determination of the mass and the moment of inertia of 
the star. 

To validate the scheme outlined above and test its ac- 
curacy, we inferred M , R and / numerically from w r and 
uj\ for three categories of compact stars, namely NS con- 
sidered in Figure CD (with AU and APRl EOS), NS con- 
sidered in lAndersson fc Kokkotasl (|1998l ) (with A, B and 
GM24 EOS), and QS (QS1 and QS2). Table[2]shows the 
percentage errors in mass, radius and moment of inertia 
of our scheme for these stars. In general, the inversion 
scheme works very well for NS and all three physical 
quantities can be determined up to a few percent or bet- 
ter. Our scheme also reproduces accurate values of M 
and / for QS. However, the inferred value of R is less 
satisfactory (~ 10 — 20%). This is understandable as the 
empirical relation only works for NS, but not QS. 
If, however, there is other independent evidence (e.g., 
via the observed thermal X-ray emission) suggesting that 
the star mig ht be a QS, we cou l d the n use equation §5§ 
proposed by iBeiger fc Haensell ()2002[ ) to infer a better 
value of R from M and /. The percentage errors SR/R 
obtained in this way are enclosed by parentheses in Ta- 
ble EU We see that the percentage error in R is decreased 
from about 10% to 1% level. In summary, we expect 
that the inversion scheme proposed in the present paper 
is applicable regardless of the nature of the compact star 
in consideration. 



5. CONCLUSION AND DISCUSSIONS 

In this paper we study the universality embedded in 
f-mode pulsations of compact stars and propose to use 
it to determine the physical parameters (including mass, 
radius and moment of inertia) of a compact star from 
which the f-mode gravitational wave signal is detected. 
We first establish a pair of empirical equations © and 
([7]) that can predict the frequency and the damping rate 
of f-modes from M and / with good accuracy for both 
NS and QS. We then apply such discovery to develop an 
inversion scheme to infer the mass, radius and moment of 
inertia of a compact star from its f-mode. In particular, 
the scheme is shown to work well for both NS and QS. 

The universality discovered in the present paper 
takes moment of inertia into the consideration, re- 
placing the role of radius used in previous attempts 
(lAndersson fc Kokkotaslll996l 119981 iBenhar et~aL1ll999l 
12004 iTsui fc Leung) l2005bD . It is worth noting that the 
radius of a star is sensitive to the low-density part of 
its EOS, whereas f-mode frequencies are expected to be 
dominated by dynamical behavior in the high-density 
regime. It is therefore physically appealing to use mo- 
ment of inertia, which measures global mass distribu- 
tion, to study f-mode oscillations. As clearly shown in 
Table [TJ the replacement of R with I indeed leads to a 
much improved universal behavior. As it is arguable that 
the accuracy of the inversion scheme directly reflects the 
quality of the universal behavior upon which the scheme 
is based, the universality discovered in the present paper 
leads to an accurate inversion scheme. 

Finally, it is worthy of mentioning that the first 
gravitational-wave search sensitive to the f-modes has 
recently been carr ied out by the LIGO detectors 
(|Abbott et al.l l2008h . While the current detectors have 
not yet detected gravitational waves directly, it is plau- 
sible that the advanced LIGO detectors, which have a 
factor of 10 improvement on the sensitivity, might de- 
tect the f-mode gravitational wave signals from compact 
stars. The inversion scheme proposed in this work will 
provide a unifying way to infer the mass, radius and mo- 
ment of inertia of both NS and QS accurately. Of course 
in reality the situation is often complicated by the pres- 
ence of magnetic field and rotation, in the future we shall 



consider the impacts of these factors on our proposal. 
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